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Abstract 

We investigate chameleonic Brans-Dicke model applied to the FRW universes. A framework to 
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that depending on the matter field and stability conditions, it is possible to realize phantom- 
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interacting with the chameleon scalar field are studied. We also compare our model with present 
day observations. 
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1. INTRODUCTION 



Recently, the observations of high redshift type la supernovae, the surveys of clusters 
of galaxies [l] — [4] , Sloan digital sky survey ( SDSS) j^J and Chandra X-ray observatory 6] 
reveal the universe accelerating expansion and that the density of matter is very much less 
than the critical density. Also the observations of Cosmic Microwave Background (CMB) 
anisotropies [?| indicate that the universe is flat and the total energy density is very close 
to the critical one The observations though determines basic cosmological parameters 
with high precisions and strongly indicates that the universe presently is dominated by a 
smoothly distributed and slowly varying dark energy (DE) component, but at the same 



time they poses a serious problem about the origin of DE 9j. A dynamical equation of 
state ( EoS) parameter that is connected directly to the evolution of the energy density in 
the universe and indirectly to the expansion of the Universe can be regarded as a suitable 
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parameter to explain the acceleration and the origin of DE [10|- 11] . In scalar-tensor theories 
interaction of the scalar field with matter ( for example in chameleon cosmology) 
24] or with geometry in Brans-Dicke (BD) cosmological models 25]— {29] can be used 
to interpret the late time acceleration. 

On the other hand, by using the well-known geometric variables, Hubble parameter and 
deceleration parameter together with the ne w g eometrical variables, the cosmological diag- 
nostic pair {r, s} ( or statefinder parameters) 30|, the acceleration expansion of the universe 
and differentiation among the cosmological models can be explained in order to better fit 
the observational data. Moreover, since the cosmic acceleration affects the expansion history 
of the universe, to understand the true nature of the driving force, mapping of the cosmic 
expansion of the universe is very crucial 31] . Hence, one requires various observational 



probes in different redshift ranges to understand the expansion history of the universe. One 
of these tests is the difference in distance modulus measurement of type la supernovae that 
helps us to testify the cosmological models. 



In here, in addition to the most of the expectations presented in 32] 33], we investigate 
the integration of both BD and Chameleon theory, which we call it " Chameleonic Brans- 
Dicke" ( CBD) model, to describe the late time acceleration of the universe and examine 
the dynamic of the universe, by applying stability analysis qualitatively. The nonminimal 
coupling term in the BD theory is equivalent to the presence of conformal transformations 
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such that one can change the form of the coupling term by applying these transforma- 
tions. Among these transformations, the ones which can transform the coupling term into 
a constant are more attractive, for they can be employed to switch between the Jordan and 
Einstein frames, which are the most discussed conformal frames 34J- 35]. Nonetheless, the 
physics is not invariant under conformal transformations, except in the weak gravitational 
field limit, and one should choose one of them as the physical frame. As the physical frame 
labels the frame we live in, it can be selected upon physical grounds consistent with prin- 
ciples and observations. In here, the CBD model in Jordan frame allows scalar field that is 
very light on cosmological scales, to couple to matter much more strongly than gravity does, 
also to couple to gravity and yet still satisfies the current experimental and observational 
constraints. The cosmological value of such a field evolves over Hubble time-scales and could 
potentially cause the late-time acceleration of our universe. From a chameleonic point of 
view, since the field mimics the background radiation/matter field, subdominant for most 
of the evolution history except at late times when it becomes dominant, it may be regarded 
as a cosmological tracker field 36[. From a BD cosmological point of view the scalar field 
is coupled to the geometry and may be regarded as an intermediate field to connect the 
field to the geometric variables in the model. In the BD model, the weak energy condition 
(WEC) imposes constraint whose validity is studied in the light of its ability to explain 
the observational data where, in par ticular, for an accelerating universe and dynamical EoS 
parameter the WEC is satisfied 37[. In CBD model where a direct coupling between the 
scalar field and the Lagrangian of the matter field is also included, the subject needs further 
studies and analysis. 

In this paper, we study the detailed evolution of the scalar field in CBD model and the 
attractor property of its solutions. Section two is devoted to a detailed formulation of the 
cosmological model. In Section three, we obtain the autonomous equations of the model and 
by using the phase plane analysis qualitatively investigate the dynamic of the system and 
the existence of a late time attractor solution. In Section four, we examine the behavior of 
the EoS parameter and deceleration parameter of the model and also perform a statefinder 
diagnostic for the model and analyze the evolving trajectories of the model in the statefinder 
parameter plane. To verify our model we also compare the distance modulus measurements 
versus redshift z derived in the model with the data obtained from the observations of type 
la supernovae. In Section five, we present summary and conclusion. 
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2. THE MODEL 



We consider the CBD gravity in the Jordan frame and in the presence of matter with the 
action given by, 

S = J d A x^-[^R - j^ a <f>, a - 2V{(j>) + 2L m f {</>)], (1) 

where R is the Ricci scalar, is the Chameleon and / or BD scalar field with the potential 
V((j)), and ui is BD parameter. We assume that f((j>) = fo(j> K and V(cj)) = Vq4> 5 where k and 
S are dimensionless constants characterizing the slope of potential V(4>) and /(0). Unlike 
the usual Einstein-Hilbert action, the matter Lagrangian L m is modified as f(4>)L m , where 
f(<f)) is an analytic function of the scalar field. The last term in the action brings about 
the nonminimal interaction between matter and the scalar field. Variation of the action (JTJ) 
with respect to the metric tensor gives, 



= ^[0,^ - ~SV0,a0 ,a ] + T^W - SW n 0] - + ^4~ T ^- ( 2 ) 

In FRW cosmology, the field equation (j2J) becomes, 

3H2= p_nj^_ 3H t + ^_m^ (3) 



On the other hand, the variation of the action (0Q) with respect to the spatially homogeneous 
scalar fields (f)(t) gives, 

9 9 3 + 2w 3 + 2u 2(3 + 2a;) ' U 

where prime means derivative with respect to <fi. We assume that the universe is filled with 
the barotropic fluid with the EoS to be p m = 7p m . From equations (J3J), ® and (jSJ), one can 
easily arrive at the modified conservation equation, 

( P J) + 3^(1 + 7 W = \0-- 3 tK/, (6) 

which readily integrates to yield, 



where M is a constant of integration. From equations (J3]) and fl4]) and in comparison with the 
standard Friedmann equations, we identify the effective energy density, p e ff, and effective 
pressure, p e ff, for the model as: 



PrnW) oTrK"<i> 2 V \<t>) 
Peff = 7 oU- + 77~l2 T~ ■ 



(8) 



Now by using equations (jHJ) and (jSJ) the conservation equation can be obtained as, 

p eff + 3H(l + u eff )p eff = 0, (10) 

with an effective EoS, p e ff = u e ffp e ff. The numerical calculation of the dynamical effective 
EoS parameter u e ff and deceleration parameter q for the model are shown in Fig 1): 



CO 





Fig. 1: The dynamical EoS parameter, u) e ff, and deceleration parameter q as a function of time. 

The functions f(<f>) = f (p K and V{<f>) = V <j) 5 , V = 20. 
For blue: « = 1, S = -4, /„ = -0.05. For black: k = 1, 5 = -2, / = 0.1 
I.Cs. 0(0) = -2, 0(0) = -1, (blue)a(O) = 2, d(0) = 0.5, (black)a(O) = 1, d(0) = 1.1 

As one can see, in Fig. 1), for k = 4, 5 = —2, the phantom crossing occurs. Also, it shows 
that the behavior of the deceleration parameter is similar to w e // ; as expected. However, 
these graphs show only the dynamical behavior of these parameters with no observationally 
validated evidence. In the following, by utilizing stability analysis of the model we revisit 
the model and further explore the observability of it by fitting the model to the experimental 
data. 
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3. PERTURBATION AND STABILITY ANALYSIS 



In this section, we study the structure of the dynamical system via phase plane analysis, 
by introducing the following dimensionless variables, 

x = ^- C = ^- v 2 = — (ID 

Using equations (jS])-©, the evolution equations of these variables become, 

, = , { z2(l±2) + , "('-y-» K + ^ + i 3(1 - 3 ^- f) (i2) 

C..«^ + ^ + |3S^^ + ^ + |^ W) 

3(l-3 7 )(2-5) 6.-12 
2(3 + 2w) X 3 + 2w ' ' 

u , 2 r 3(l -3 7 )(2-<5) 3(1 + 7). 2 35-6 2l 
" = < + 2 C + 1 2(3 + L) + - V 2 '^ + 3T2^" >' (14) 

where prime in here and from now on is taken to be derivative with respect to N — Ina , 
together with the Friedmann constraint equation d3J) which becomes 

x 2 -C + ^C 2 -v 2 = i- (15) 

In term of the new dynamical variable, we also have, 

IL = o t - V _ r 3(l-37)(2-J) _ 3(l+ 7 ) 2 _ 38-6 2 

H 2 ^2^ [ 2(3 + 2w) 2 lX 3 + 2u ' ' 1 J 

Using the constraint (fl5l) . the three equations ( fT2l - (fT4l) now reduce to the following two 
equations: 

* - *{^> + l^f^K + + !fc^> (17, 

35-6, , 2 l. w a2n , 3(l-3 7 )(2-5) 2 65-12. 1 2 . u; 9 , 

It is more convenient to investigate the properties of the dynamical system equations (ITT)) 
and (ITS)) than equations f|T2|) - fri4|) . by calculating the critical points ( or fixed points) we 

6 



study the stability of these points. Critical points are always exact constant solutions in the 
context of autonomous dynamical systems. These points are often the extreme points of 
the orbits and therefore describe the asymptotic behavior of the system. In the following, 
we find fixed points by simultaneously solving x' — an d C = 0- Substituting linear 
perturbations x' ~ > x' + $x'i C ~ * C + about the critical points into the two independent 
equations (ITT)) and ( Tl8|) . to the first orders in the perturbations, yields two eigenvalues 
= 1,2). Stability requires the real part of the eigenvalues to be negative. Solving the 
above equations we find seven fixed points which some of them explicitly depend on 7, k 
and S, as illustrated in tables I to III, for 7 = 0, 7 = 1/3 and 7 = — 1 respectively. 



TABLE I: critical points for 7 = 



points PI P2 


P3 


P4 


P5 


P6 


P7 


x 





A 


A 


B 


B 


V3(k+8) 


V3{k+8) 


-K+4S 


-K+45 


C 3 + VIE 3- VIE 


4-25 


2k-A 


2k-4 


12 


12 


5+3 


K + 8 


K+8 


-K+45 


-K+45 



where A = V~20k - 280 + k 2 and B = V-165 2 + 565 + ASk - 8k + 72 



TABLE II: critical points for 7=3 



points PI P2 P3 P4 P5 


P6 


P7 








x 0-11 


V5' 2 -45~ 
5 


6 V<5 2 -4<5-6 
.5 


C 3 + VIE 3 VIE 4 5 ~ f 


4 
5 


4 
5 


TABLE III: critical points 


for 7 = 


-1 


points PI P2 P3 P4 


P5 


P6 P7 


x 


1 5-2 


c c 

3k+9 3k+9 


C 3 + VIE 3 VIE 4 5 " 25 





4-2/t 4-2k 

K+3 K+3 
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where C = V165 - 60k - 15k 2 . 



In the following we will investigate the stability of the model with respect to the above 
three specific choices of the EoS parameter for the matter in the universe, i.e. 7 = 0, 7 = 1/3 
and —1. 



3.1. Stability for 7 = 

In the case of 7 = 0, the stability properties of the seven critical points in this system 
are shown in the following: 

PI : Ax = (3 + v / l5)(-2 + Vl5 + 5), A 2 = ^(3 + v / 15)(-5 + + re), stable for 

8 

5 < 2 - re < 5 - Vl5 

P2 : Ai = (3- v / 15)(-2- VTE + 6), A 2 = ^(3 - v / 15)(-5 - Vl5 + re), stable for 

8 

5 > 2 + >/l5, re > 5 + Vl5 
P3 : Ai = (-11 - 45 + 5 2 )/(3 + 5), A 2 = (45 2 - 145 - 5re + 2re - 18)/(12 + 45), stable for 
i)5 < -3, s < k , m)2 < 5 < 2 + \/l5, £ < k m)2 - < <5 < 2, re < £ UH)Vk, 5 = 2 
PA, P5 : Ai = O 2 - Ak - 56)/(4re + 32), A 2 = (k 2 + Ak - A5k + 85 + 48)/(2re + 16), stable for 
i)re < -8, pt< 5 , ii)2 < re < 2 + \/60, /i < 5 m)2 - < re < 2, 5 < /1 



P6, P7 : Ai 2 = — — — stable for some values of 5 and re, (19) 

—10k + 40o 

where e = 2(25 2 - 75 - 9)/ (5 - 2), /1 = (k 2 + 4k + 48)/(4k - 8), and 

Q = -15005 2 + 226805 - 420k5 + 26820 + 180k + 525k 2 + 7805k 2 + 30k 3 5 - 60k 3 + 4805 3 k - 

9605 3 - 1920k5 2 - 2405 2 k 2 . 

We see that all the critical points are stable for the given conditions on the stability 
parameters re and 5. 

In Fig. 2), the attraction of trajectories to the critical points PI to P7 in the phase plane 
is shown for the given conditions on re and 5. Fig. 2a), belongs to the critical point PI. 
Fig.2b) is an example of the stability of the critical point P3 that satisfies the condition i). 
In Fig 2c) since there are two critical points P4 and P5 with the same stability conditions, 
in here, condition i) gives us two stable critical points. The same argument is applies to the 
Fig 2d) with respect to the points P6 and P7 where for some numerical values of 5 and re 
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the two eigenvalues become negative and critical points become stables. In addition, since 
for the critical point PI and P3, x = an d C 7^ for the stability conditions on 5 and k, 
these points corresponds to a scalar field dominated era. 



6=-2,k=-6 






X 



8=10,k=-10,h= 



26 
3 





5=-2,k=3 



4- 

C 3 " 




unstable,^ J 



Fig. 2: The attractor property of the dynamical system in the phase plane for the case ' 
The plots a) to d) are respectively examples of the critical points PI to P7. 
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3.2. Stability for 7 = ± 



In the case of 7 = |, the stability properties of the seven critical points in this system 
are shown in the following: 



5 \/15 

P1,P2 : Ai = (3± v / 15)(2±v / 15 + 5), A 2 = - ± unstable 

P3 : X 1 = (-11 -45 + 5 2 )/{3 + 5), A 2 = {5 2 -45- 6)/(3 + 5), stable for: 
z)Vk, 5 < -3 , zi)V«, 2 - VTo < 5 < 2 + vTo 



(20) 



P4,P5 



Ai = 4, A 2 



2 = — 1, unstable (saddle point) 



P6,P7 



2 - <5 ± V100 + 605 - 155 2 
25 



, stable for z')Vk, 5 > 2 + v 7 !!), zz)Vk, 5 < 2 - \/l0 



We explicitly find that the critical points P3, P6, P7 are stable for the given conditions on 
the stability parameters k and 5. In Fig. 3), the attraction of trajectories to the critical 
points in the phase plane is shown for the given conditions on k and 5. In Fig. 3a), although 
the given stability parameter 5 = 4 satisfies the condition i) of the critical point P3, since 
it also satisfies the conditions ii) of the critical points P6 and P7, as can be seen, there are 
three stable points displayed in the graph. In addition, since for the critical point P3, x = 
and ( 7^ for all 5 and k, this point corresponds to a scalar field dominated era. 
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k=1,8=-4 



k= 10,8=0 




k=1,8=10 



C 3 
2 
1 ■ 



\ \ \ \ \ \ 
' , \ \ \ \ 




stable • 

. unstable* 




stable 
unstable* 



s/ J J- ^ \ \ \ \ % \t 

^ W V 1 1 ^ \ \ "\ "\ 



o 

X 

k=10,8=-2 




/ / / / 

stable • 
unstable* 



3: The attractor property of the dynamical system in the phase plane for the case 7 
The plots a) to d) are respectively examples of the critical points P3, P6 and P7. 
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3.3. Stability for 7 = 



-1 



In this case, we find the stability conditions for the seven fixed points in the system as 
shown in the following: 

PI : Ai = (3 + v / 15)(2 + VIE + 5), A 2 = (3 + v / 15)(2 + VIE+k), stable for 

8 < 2 - VIE, k<2- VIE 
P2 : Ai = (3- v / 15)(2- VIE + 5), A 2 = (3 - v / 15)(2 -VIE+k), stable for 

5 < 2 - VIE, k<2- VvE 

(-11-45 + 5 2 ) {-25 + 2n + 5 2 - k5) 
P3 : A, = _ , A 2 = ^ , stable for 

i)8 < k, 8 < -3 ii)8 < n, 2 < 5 <2 + VIE Hi) 8 > k, 2 - VIE < 5 < 2 
3 1 

P4, P5 : Ai >2 = ~ ± — V705 - 240k + 1205k - 2405, stable for i)8 < 0, re > a ii)8 > 0, re < a 

iii)8 < 2, k > (3 iiii)8 > 2, k < f3 

-4k + k 2 -11 , -2(-k 2 + 2k + k5 - 28) ln „ 
P6, P7 : Ai = , A 2 = i 1 L : stable for 

K + 3 K + 3 

z')k < 8, 2 < k < 2 + VIE U)k > 8,2- V^lo" < k < 2 

where alpha = 2 ^ + ^ , (3 = |^-^z^^ 

Again, we explicitly find that all of the critical points are stable by the given conditions 
on k and 8. In Fig. 4), the attraction of trajectories to the critical points PI and P3 
in the phase plane is shown for particular choice of k and 8. In Fig. 4b), similar to the 
argument we had in case 7 = 1/3, although the given stability parameters 8 = —4, k = 10 
satisfies the condition i) of critical point P3, since it also satisfies the conditions iii) of the 
critical points P4 and P5, there are three stable points displayed in the graph. Also, Fig. 
5) shows the attraction of trajectories to the critical points P4 to P7 in the phase plane 
for the conditions on k and 8. Again, the extra stable points shown in Fig. 5 a), c), d), 
and f) correspond to the points satisfied by the stability conditions for other critical points. 
Moreover, since for the critical points PI, P2 and P3, \ = an d C 7^ for all 8 and k, these 
points correspond to a scalar field dominated era. Similarly, for the critical points P4 and 
P5, C = and x 7^ for all 8 and k, therefore, these points correspond to either a matter 
dominated era or /(</>) dominated era. 
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unstable • / / i j i i \ \ \ \S\\s, 



stable 



-l 
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k =6, S =5 
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\ \ / / y ^».-^- ...... 
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\ \ 
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unstable • 



-2 



' ""^ft \ \ \ I I 




-10 ,8 =0 



d) f f 

\ t f f 
t t t f 
t f / / 



\ \ \ \ t 
\ \ \ \ \ 
WW \ 
WW \ 
w \ \ \ 
W \ \ \ 
\ X t f 
\\\ t f f 
\ \ t t t 




stable 
unstable 



Fig. 4: The attractor property of the dynamical system in the phase plane for the case 7 
The plots a) to d) are respectively examples of the critical points PI, P2 and P3. 
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k =5,5=-2,orf 



7- 
6- 
5- 
4 
3- 
2- 
1- 
0- 
-1- 
-2- 



J J S if 

\ I I / / 
\ l t / /rf? 
\ \ I / 



\ *\ \ \ \ \ 

\ \ \ \ ^ 

"N, \ \ \ \ | 

s> \ \ \ i i 

k\ \ \ \ I 

* 1 i 




\ \ \ 

unstable • \ \ [ 



stable 



k —5,5—2, oc. — 1 

f f f / / / y / \\ \ \ \ ^ ♦ f f f 
, f t t / .-»— -"x \ / V \ \ \ T f f 

\ t f / /J^^Kv \ \ j / Z^^nK \ t f f 

t t f I ) \ \ t 

\ t t 

\ t '// \. \ W / ^ - - \V f / 

" t t 

\ \ f/^\.\l 7 ^ Ml / / 

\ \ 

/ / 
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unstable • 




-l 



K=10,8 =0,(3= 



47 
16 



K=0,6=3,P= 



WW \ \ 

\ \ \ \ 

N» \ \ \ \ 

\ \ \ 




stable 
unstable 




unstable 



k=5,6 =6 



K=0,S =-6 




f) 1 1 1 J / / / /i 

..III . 1 



)\\\ \ ^ ^ | ^ 
\ \ V 




stable 
unstable 



5: The attractor property of the dynamical system in the phase plane for the case 7 = 
The plots a) to f) are respectively examples of the critical points P4, P5, P6 and PI. 
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4. COSMOLOGICAL PARAMETERS: EOS PARAMETER, STATEFINDERS 
AND DISTANCE MODULUS 



In order to understand the behavior of the universe and its dynamics we need to study 
the cosmological parameters such as EoS parameter, deceleration parameter and statefmder. 
We also need to verify our model with the current observational data which in here we test 
the distance modulus. Some of these parameters like EoS and deceleration parameters or 
even distance modulus analytically and/or numerically have been investigated by many 
authors for variety of cosmological models. Utilizing stability analysis enables us to have 
a deeper understanding of the problem and best fit the model with the observational data. 
For our model the effective EoS and deceleration parameters are defined by q = — 1 — M% 
and u e fj = — 1 — where in terms of new dynamical variables is given by equation 
(TIB"]) . The statefmder parameters are another set of parameters that are studied in here. 
These parameters are defined by r = — 3q — 2 and s = 3 ^~ 1 1 ^ where q already defined in 
terms of new dynamical variables and we have = ^(jJt) + 2(-j|) 2 . 

In tables IV to VI for three different scenarios, 7 = 0, 1/3, —1, the above cosmological 
parameters are calculated as: 



TABLE IV: Properties of the critical points for 7 = 
points q uj e fj r 



PI 


5 + ^ 


3+ fv^ 


^ + fv / 15 


P2 


5 - 


3- §\/l5 


W -f v^5 


P3 


S 2 -4S-1 


2<5 2 -9<5-5 


B 


3+5 


9+35 


(6+3)* 


P4,5 


k 2 -4k+24 


k 2 -6k+8 


k 4 -8k 3 +76k 2 -320k+ 


32+4k 


48+6k 


8(k+8) 2 


P6,7 


k 2 -36<5k+56<5 2 -16k- 176<5-72-3<5k 2 +12<5 2 k 


^6,7 




5(-k+4<5) 2 



17V15+67 



2V15+9 6+V15 
17^-67 1 
2VTE-9 6-V15 

<5 2 -35+2 



c 



<5+3 



3 -6k 2 +56k-336 32+4k 
6k 2 +24k-192 k 2 +56 



«6,7 



1 

9+1 



where 

-3k 2 -325k+325 2 -32k-3525-144-6<5k 2 +245 2 k) 
W 6,7 - 15(-k+45) 2 , 

_ (2160-81445 2 + 11925k-5765 3 +76k 2 +4305k 2 -24805 2 k-11k 3 +3(5k 3 - 204.5 2 k 2 +768.5 3 k-168k+16512(5) 
r 6J — 5(_ K+ 45)3 

_ 2(2160-81445 2 +11925k-8965 3 +76k 2 +370(5b 2 -22405 2 k-6k 3 +3(5k 3 -2045 2 k 2 +7685 3 k-168k+165125) 
S 6> 7 ~~ 3(-k+4<5)(-3k 2 -32<5k+32<5 2 -32k-3525-144-6<5k 2 +24<5 2 k) ' 

where D = 9 - 3645 + 105 5 + 3155 2 + 565 3 - 745\ and E = gi^zg^ggt^^ , 
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TABLE V: Properties of the critical points for 7=3 



points q w e// 



PI 5 + 3 + %y/lE ^3 + 51 y^5 



17V15+67 



P2 5-VT5 3-1V15 ^-4f^l5 



2 ' 2 v 2^+9 6+^ 

2./TE 203 51. /Te 17^-67 1 



3 v 2 2 v 2 v / 15-9 6-v^5 



po 5 2 -45-l 25 2 -95-5 D p 5 2 -35+2 

- r ° 3+5 9+35 (5+3) 3 ^ 5+3 

P4.5 1 



1 o 4 1 

3 ° 3 2 



Pfi 7 ^2 ^4 (-465+76+35 2 ) 4 (~235+38+5 2 ) 5 

rD >' 5 35 P 3 5(5-4) 25-2 



TABLE VI: Properties of the critical points for 7 = — 1 



points 


q 








r 




s 




P 


PI 






3 + §\/l5 


203 


+ ^\ 


15 


17^4 


-67 


1 




2 


2^4 


-9 


6+^ 


P2 


5 - y/TE 




3- |Vl5 


203 




15 


17\/i5- 


-67 


1 




2 




2\/l5- 


-9 




P3 


5 2 -45-l 




25 2 -95-5 




D 








3+5 


3+5 




9+35 




(5+3) 3 






5 2 -35+2 


P4,5 


-1 




-1 




1 









OO 


P6,7 


(-7k 2 +13k+3-k 3 +8k5) 


K- 


13k 2 +32k+15-2k 3 - 


hl6fc5) 


F 




G 




-(k+3) 2 


(k+3) 2 


3 


(«+3)« 






(- 
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where F — ~(~ 9 ~ 112k2 ^~ 6k4 + 264k ^~ 160k ^ 2 + 80k2 ^ 2 + 32k3 ~ 2k5 + 12k ~ 99k2 ) 

(k+3) 3 ' 
, „ _ |(18-112k 2 5-6k 4 +264k5-160k5 2 +80k 2 5 2 +33k 3 -2k 5 +39k-90k 2 ) 
and G — (-13k 2 +32k+15-2k 3 +16k5)(k+3) ' 

In the above tables p is the power of time t in a(t) oc t p obtained in terms of the dynamical 
variables. For 7 = 0, from tables I and IV, the critical points PI and P2 are independent 
of the stability parameters k and S, having numerical values. These critical points show a 
decelerating universe with positive effective EoS parameter. On the other hand, the critical 
point P3 is stable under some conditions on 5 while independent of k while the points P4 
to P7 are both 5 and k dependent. All these critical points can be regarded as attractors 
for all the trajectories in the statefinder phase plane if satisfy the stability condition. In 
case of 7 = 1/3, as can be seen from tables II and V, the unstable critical points PI, P2, 
P4 and P5 are independent of the stability parameters k and 5, given the numerical values. 
On the other hand, the critical points P3 , P6 and P7 are stable under some conditions on 
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S while independent of k and so can be regarded as attractors for all the trajectories in the 
statefinder phase plane. For 7 = — 1, from tables III and VI, while the critical points PI, 
P2, P4 and P5 are dependent of the conditions on k and 5, their physical parameters have 
numerical values. On the other hand the critical points P3 , P6 and P7 are independent of 
k. All these points can be regarded as attractors for the trajectories in the statefinder phase 
plane. 

We can discuss the dynamics of the EoS parameter and deceleration parameter in Fig. 
6) in the case of 7 = 0. As is shown, for a variety of stability parameters, the evolving 
effective EoS parameter, w e // ; and the deceleration parameter q are numerically computed. 
From the graph, the EoS parameters for the stability parameters 5 = 6, k = 3 (green curve) 
do not hit the phantom divide line , while for 8 — 7, k — 1 and 5 = 6,k = 3 ( red and 
blue curves) the trajectories cross the line. It also shows a symmetrical behavior in the 
parameters about ln{a) = 0.15 ( z = —0.139). Moreover, the deceleration parameter, q at 
the present time is negative and between about —0.9 and —2.1 (an accelerating universe) 
for the given stability parameters. For both blue and red curves, the crossing occurs at 
about ln(a) = —0.2 ( z = 0.221) in the past while in the future at ln{a) = 0.35 (z = —0.3) 
for blue curve and ln(a) = 0.5 (z = —0.393) for red curve. Using the above argument, and 
from the Fig. 13)left) one can see that all the trajectories crossing phantom line in the 
interval z = — 1 are best fitted with the observational data while the red curve in the 
interval z = 1 — 2 better fits the data. Noting that all the trajectories satisfy the stable 
critical points P6,7 condition. With the selected stability parameters in blue trajectory, 
uj e ff = —1/9 and in the red and green trajectories (J e // = —1/21 which in both cases it 
occurs sometimes in the far past and future, as can be seen in the graph. 
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Fig. 6: The evolution of the effective EoS parameter u e ff and deceleration parameter 
as a function of z for 7 = 0. Stability conditions: (blue)<5 = 7,k = l-(green & red)<5 = 6,k 
I.Cs.:(blue & green)x(O) = 0.1, C(0) = -0.5 (red)x(O) = 1/3, C(0) = -2/9 

In Fig. 7), the evolving effective EoS parameter, u) e ff, and the deceleration parameter q 
are shown for 7 = 1/3. It can be seen that the EoS parameter for the stability parameter 
5 = 4 (blue curve) does not hit the phantom divide line while for 5 = 6 (black curve) crosses 
the line. It also shows a symmetrical behavior of the curves sometimes in the future about 
lines ln(a) = 0.2 ( z = —0.18) and ln(a) = 1.5 ( z = —0.78) respectively for the above 
two stability parameters. Moreover, the trajectory in the case 5 = 4 while symmetrically 
and sinusoidally oscillates, eventually diminishes in far past and future. The graph also 
shows that the deceleration parameter, q, which has a similar behavior to the effective EoS 
parameter, is at the present time negative and —0.6 or —1.4 (an accelerating universe) for 
the given stability parameters. With respect to these two different EoS trajectories and 
stability analysis, a question arises as to the validity of the solutions since one of them is 
crossing the phantom line and the other one not. The answer is that the trajectory which 
best fit the observational data should be considered as physical one. From Fig. 13)middle), 
as will be discussed, one can see that the difference in the distance modulus is best fitted 
with the observational data for the trajectory that crosses the phantom line. Noting that 
both trajectories correspond to the stable critical points, the blue one satisfies the critical 
points P3)ii) condition, while the black one satisfies the P6,7)i) condition. With the selected 
stability parameter in blue trajectory, u e ff = —0.5 and in the black trajectory u e ff = 1/9 
which in both cases it occurs sometimes in the past and future as can be seen in the graph. 
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Fig. 7: The evolution of the effective EoS parameter u e ff and deceleration parameter q 
as a function of ln(a) for 7 = 1/3 (blue):5 = 4, (black):5=6 - I.C. x(0) = 0.21, C(0) = -0.4 

With the same argument as above, we can discuss the dynamics of the EoS parameter 
and deceleration parameter in Fig. 8) in the case of 7 = — 1. As is shown, for a variety 
of stability parameters, the evolving effective EoS parameter, oo e ff, and the deceleration 
parameter q are numerically computed. From the graph, the EoS parameters for the 
stability parameters 5 = 6, k = 3 do not hit the phantom divide line (green and black 
curves), while for 8 — 7, k — 1 and 5 = 5, ft = 5 the trajectories cross the line ( red and 
blue curves). It also shows a symmetrical behavior in the parameters about ln(a) = 0.1 ( 
z = — 0.1)and ln(a) = 0.2 ( z = —0.18). Moreover, the deceleration parameter, q at the 
present time is negative and between about —1 and —2.3 (an accelerating universe) for 
the given stability parameters. For the stability parameter 8 — 7, k = 1 (blue curve), the 
crossing occurs at ln(a) = 0.3 ( z = —0.26) in the future and ln(a) = —0.1 ( z = 0.11) in 
the past. For the stability parameter 5 = 5, k = 5 (red curve), it occurs at ln(a) = 1.2 ( 
z = —0.7) in the future and ln(a) = — 1.2 ( z — 2.3) in the past and then again become 
tangent to the line in the future and past. Using the above argument, and from the Fig. 
13)right) one can see the blue trajectory crossing phantom line is best fitted with the 
observational data in comparison with the other trajectories. Noting that the red trajectory 
correspond to a critical point that is not stable, the green and black trajectories both 
with the same stability parameter and different initial conditions satisfy the critical points 
P6,7)i) condition, while the blue trajectory satisfies the P3)iii) condition. 
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ln(a) ln(a) 

Fig. 8: The evolution of the effective EoS parameter ui e ff and deceleration parameter q 

as a function of ln{a) for 7 = — 1. Stability conditions: (blue) 5 = 7,k = 1 
-(red)<5 = 5,k = 5-(green & black)<5 = 6,« = 3 I.Cs.:(blue & green)x(O) = 0.1, C(0) = -0.5 
(red)x(O) = 1/3, C(0) = -2/9-(black) X (0) = 3/10, C(0) = -1/3 

As a complementarity, Fig. 9) shows statefinder diagrams for {s, q} and {r, q} evolu- 
tionary trajectories for 7 = 0. From this figure, we clearly see the position of the current 
values of the statefinder parameters for the two given stability conditions with respect to 
the LCDM, SCDM and SS states. The trajectories commence evolving from the same state 
and end their evolution to the same common state. 
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Fig. 9: Trajectories in the statefinder plane {r, q} and {s, q} for the 
model with 7 = O.The arrows show the direction of the time evolution. 
Stability conditions: (blue) S = 7,k = l-(green & red)5 = 6,k = 3 
I.Cs. (blue & green)x(O) = 0.1, C(0) = -0.5 (red)x(O) = 1/3, C(0) = -2/9 

Also, Fig. 10) shows statefinder diagrams for {s, q} and {r, q} evolutionary trajectories 
for 7 = 1/3. From this figure, we clearly see the position of the current values of the 
statefinder parameters for the two given stability conditions with respect to the LCDM, 
SCDM and SS states. The graph also shows that in both cases the trajectories evolve to 
an unstable critical point. They commence evolving from the same state and end their 
evolution to the same common state. 
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Fig. 10: Trajectories in the statefinder plane {r, q} and {s, q} for the 
model with 7 = 1/3. The arrows show the direction of the time evolution. 
Stability conditions: (blue) S = 4-(black)<5 = 6. I.Cs. x(0) = 0.21, C(0) = -0.4 

in the case of 7 = — 1, Fig. 11) shows statefinder diagrams for {s, q} and {r, q} evolu- 
tionary trajectories. From this figure, the position of the current values of the statefinder 
parameters for the four given stability conditions with respect to the LCDM, SCDM and 
SS states are shown. Likewise the other cases, the trajectories commence evolving from the 
same state in the past and end their evolution to the same common state. 
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Fig. 11: Trajectories in the statefinder plane {r, q} and {s, q} for the model with 7 = — 1. 
The arrows show the direction of the time evolution. Stability conditions: (blue) 5 = 7,k = 1 
(green)5 = 6,k = 3-(red)5 = 5,k = 5-(black)<5 = 6,k = 3 I.Cs. (blue & green)x(O) = 0.1, £(0) = 
(red)x(O) = 1/3, C(0) = -2/9-(black) X (0) = 3/10, C(0) = -1/3 

In Fig. 12), the statefinder diagrams {r, s} for 7 = 0(left), 7 = l/3(middle) and 
7 = —1 (right) are shown. Again the current state of the universe for the given stability 
conditions with respect to the LCDM scenario is presented. The graph shows that in case of 
7 = 1/3 the trajectories start from an unstable critical point and end to the same unstab 
point which may be indirect evidence of a cyclic universe predicted by string theory 



The current values of the trajectories in our model for 7 = 0,1/3,-1 and their distance 
from LCDM scenario can be seen explicitly in these diagrams. Note that the trajectories of 
our model will pass through LCDM fixed point. It is of interest to see that the trajectories 
have different forms before reaching the LCDM state 30] . This behavior demonstrated that 
the statefinder can successfully characterize and differentiate between various DE models. 
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Fig. 12: Trajectories in the statefinder plane {r, s} for 7 = O(left), 7 = ^(middle), 
7 = — l(right). For 7 = 0:(blue)5 = 7,k = l-(green &red)<5 = 6, ft = 3 
I.Cs.:(blue & green)x(O) = 0.1, C(0) = -0.5 (red)x(O) = 1/3, C(0) = -2/9. 

For 7 = §: (blue)5 = 4-(black)<5 = 6. I.Cs. x(0) = 0.21, C(0) = -0.4 
For 7 = —1: (blue)5 = 7,k = l-(red)5 = 5,k = 5-(green & black)5 = 6,k = 3 
I.Cs. (blue & green)x(O) = 0.1, C(0) = -0.5 (red)x(O) = 1/3, C(0) = -2/9 
-(black)x(O) = 3/10, C(0) = -1/3 

Finally, to validate the model, we check it against observational data from our past light 
cone such as distance modulus which is defined by fi(z) = m(z)—M, where m is the apparent 
magnitude of the source and M its absolute magnitude. The distance modulus is related to 
the luminosity distance via fi(z) = 5log w d,L(z) + 25 where d^z) is measured in units of 
megaparsecs. In Fig 13), the distance modulus obtained from observation at high redshift 
is compared with the one numerically computed by means of the model using the stability 
analysis results. From Fig. 13)left) for 7 = 0, as can be seen the n(z) for three sets of 
stability parameters 5 and 5 in the model are compared with the observational data. The 
curves are in general in the best fit with the observational data in the interval z = — 1 while 
for the red curve this continues to the point z = 2. From Fig. 13)middle) for 7 = 1/3, as 
can be seen the /i(z) for two stability parameters 5 = 4 and 5 = 6 in the model are compared 
with the observational data. In 5 = 4 case (blue curve) where the critical point p3 is stable 
(with q = —1/7 and u e ff = —9/21), the graph shows slightly better fit the observational 
data in comparison with the case 5 = 6 (black curve) where the model is unstable. 

In Fig.l3)right) for 7 = —1, as can be seen, the fi(z) for four different stability 
parameters are compared with the data. For the stability parameter 5 = 5, k = 5 (red 
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curve) which corresponds to an unstable universe, the numerically calculated fi(z) is out of 
range. For the stability parameter * = 6, k = 3 (green curve) where the model contains 
stable critical points at some stage and the universe passes through the LCDM state for 
the given initial conditions the model fit the observational data. The black curve has the 
same stability parameters but with the given initial conditions it does not fit the data. For 
the stability parameter * = 7, k = 1 (blue curve) the model contains stable critical points 
and the universe at some stage passes through LCDM state and as shown the model best 
fit the observational data. 




Fig. 13:The graph of the distance modulus fj,(zO for 7 = O(left), 7 = | (middle), 7 = — l(right). 
For 7 = 0:(blue)* = 7,k = l-(green &red)* = 6,k = 3 
LCs.:(blue & green)x(O) = 0.1, C(0) = -0.5 (red)%(0) = 1/3, C(0) = -2/9. 

For 7 = §: (blue)* = 4-(black)* = 6. I.Cs. x(0) = 0.21, ((0) = -0.4 
For 7 = —1: (blue)* = 7,k = l-(red)* = 5,k = 5-(green & black)* = 6,k = 3 
I.Cs.(blue & green)x(O) = 0.1, C(0) = -0.5 (red)x(O) = 1/3, C(0) = -2/9 
-(black)x(O) = 3/10, C(0) = -1/3 



5. SUMMERY AND CONCLUSION 

This paper is partially designed to study the attractor solutions of the CBD cosmology 
by utilizing the 2-dimensional phase space of the theory. The potential V(4>) and coupling 
scalar function /(</>) in the model is considered to be in power law forms in the phase space. 
The matter Lagrangian in the model is regarded as a perfect fluid with three kinds of EoS 
parameters, 7 = 0,1/3,-1. The stability analysis gives the corresponding conditions for 
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tracking attractor and determines the universe behavior in the past and future. By stability 
analysis, the critical points in the model for the matter fields are evaluated to be stable if 
satisfying stability conditions, Figs. 2-5. 

We then study the cosmological parameters such as effective EoS parameter, oo e ff, decel- 
eration parameter, q, and statefinder parameters for the model in terms of the dynamical 
variables introduced in the stability section. The results show that depending on the stabil- 
ity parameters, the model produces the expansion of the universe characteristic of phantom 
cosmology, despite the absence of phantom energy. It also shows that the EoS parameter 
may cross the cosmological divide line depending on the stability conditions in the past and 
future in all 7 = 0, 1/3, —1 cases. Then we perform a statefinder diagnostic for the CBD 
model. It is shown that the evolving trajectories of different scenarios in the {r, s} plane 
is different from those of other DE models. From the graphs we also find that in case of 
7 = 1/3 the universe may have a cyclic behavior between two unstable states, while passing 
through LCDM state on its way. 

In conclusion, by fitting the observational data to the model, we compare the distance 
modulus-redshift relation for 7 = 0, 1/3, —1 cases, depending on the stability conditions. In 
general, in comparison of the observational data with the model for different 7's and stability 
parameters show that the scenarios with phantom crossing better fit the data. Among those 
scenarios with phantom crossing and best fit with the observational data, the 7 = case for 
the stability parameters 5 = 6 and k = 3 ( red curve) is more interesting because the from 
the Fig. 6) we see that the effective EoS parameter, co e ff, and deceleration parameter q start 
from dust-like dominated era and followed by the dark energy dominated era at the present 
time. The symmetric behaviour of the graphs shows that the universe continues being in 
the energy dominated era for a while and finally reaches the dust-like dominated era again. 
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